% 子函数：用文档1.2.2节近似公式计算相位（φₘ(t)≈2πft + εₘ + A_Rcos(2π(f_Rt+m/M)+α_R)）
% -------------------------------------------------------------------------
function approxPhase = calcApproxPhase(timeSeg, psiMSeg, rotationFreqSeg, config)
    % 输入：
    %   timeSeg - 单段时间向量(s)
    %   psiMSeg - 第m天线的旋转相位ψₘ(t)(rad)
    %   rotationFreqSeg - 实时旋转频率f_R(t)(Hz)
    %   config - 信号配置（含A_R、α_R、εₘ等参数）
    % 输出：
    %   approxPhase - 近似相位(rad)
    
    % 1. 提取近似公式参数（从config读取或设默认值）
    f_common = config.signal.CenterFrequency; % 公共载波频率f（文档1.2.2节定义）
    eps_m = config.antennas.phaseOffsets(1); % 第1天线独立相位εₘ(rad)
    A_R = config.rotation.amplitude; % 旋转多普勒幅度A_R(rad)（config中已定义）
    alpha_R = config.rotation.initialPhase; % 旋转相位α_R(rad)（config中已定义）
    m = 1; % 当前天线索引（第1天线）
    M = config.antennas.count; % 总天线数M
    
    % 2. 文档1.2.2节近似公式：φₘ(t)≈2πft + εₘ + A_Rcos(2π(f_Rt + m/M) + α_R)
    approxPhase = 2 * pi * f_common * timeSeg + eps_m + ...
                  A_R * cos(2 * pi * (rotationFreqSeg .* timeSeg + m/M) + alpha_R);
    
    % 3. 相位转换到-pi, pi
    approxPhase = wrapToPi(approxPhase);
end